Oxygen-resistant hydrogenases and methods for designing and making same

ABSTRACT

The invention provides oxygen-resistant iron-hydrogenases ([Fe]-hydrogenases) for use in the production of H 2 . Methods used in the design and engineering of the oxygen-resistant [Fe]-hydrogenases are disclosed, as are the methods of transforming and culturing appropriate host cells with the oxygen-resistant [Fe]-hyrdogenases. Finally, the invention provides methods for utilizing the transformed, oxygen insensitive, host cells in the bulk production of H 2  in a light catalyzed reaction having water as the reactant.

GOVERNMENT INTEREST

The United States Government has rights in this invention pursuant to contract DE-AC36-99G010337 between the United States Department of Energy and the National Renewable Energy Laboratory, a division of Midwest Research Institute.

FIELD OF THE INVENTION

The present invention relates to hydrogen-production by microorganisms, for example hydrogen production by green algae. More specifically, the invention relates to methods for designing and engineering hydrogenase enzymes with improved oxygen resistance, and to the methods for transforming microorganisms to express these oxygen-resistant hydrogenase enzymes for use in the production of hydrogen in) an oxygen containing environment.

BACKGROUND OF THE INVENTION

Hydrogen (H2) is becoming an attractive alternative energy source to fossil fuels m due to its clean emissions and potential for cost effective production by microorganisms. As such, microorganisms that metabolize H₂ are being investigated for their potential use in H₂-production. A microorganism of particular interest for H₂ production is the green alga, Chlaydomonas reinhardtii, which is able to catalyze light-dependent, H₂ production utilizing water as a reductant. Ghiradi et al., (2000) Trends Biotech. 18(12):506-511; Melis et al., (2001) Plant Physiol. 127:740-748. The benefits of using an algal system for H₂-production include the use of renewable substrates (light and water) and its potential cost-effectiveness. Melis A, Int. J. Hyd. Energy 27:1217-1228. As such, there is a great deal of interest in optimizing H₂-production by green algae to maximize the potential benefit as an alternative energy source.

Chlaydomonas reinhardtii, and other like microorganisms, are able to express a class of H₂ metabolizing enzymes called hydogenases. Members of this enzyme family function in either H₂-uptake (as a means to provide reductant for substrate oxidation) or H₂-production (as a means to eliminate excess reducing equivalents). Characterization of various hydrogenases from multiple organisms has identified three principle hydrogenase types, broadly classified by the chemical nature of their active sites: [Fe]-hydrogenase, [NiFe]-hydrogenase, and non-metallic (organic) hydrogenase. Vignais et al, (2001) FEMS Micro. Rev. 25:455-501; Adams M. W., Biochem. Biophys. Acta. 1020:115-145; Buurman et al., (2000) FEBS Letts. 485:200-204. More particularly, [Fe]-hyrdogenase have an active site containing a [4Fe-4S]-center bridged to a [2Fe-2S]-center (H-cluster) (Peters et al., (1998) Science 282:1853-1858; Nicolet et al., (1998) Structure 7:13-23), and the [NiFe]-hydrogenase have an active site containing a [4Fe4S]-center bridged to a [NiFe]-center (Volbeda et al., (1995) Nature 373:580-587). Coordination of the metal prosthetic groups to the active sites is made by cysteinyl, CN⁻, and CO ligands. Further, within each hydrogenase group are monomeric, or multimeric enzymes, that can be either cytoplasmic or membrane bound within the cell. Vignais et al., Supra.

Although there are differences within the active sites between different families of hydrogenase, as well as between the subunit composition and localization between hydrogenase families, most, if not all studied hydrogenases have exhibited some degree of sensitivity to inhibition by CO and O₂. Adams M. W. W; Volbeda et al., (1990) Int. J. Hyd. Energy 27:1449-1461. Hydrogenase sensitivity to these inhibitors correlates to some degree to the type of prosthetic group that forms the active site, for example, [Fe]-hydrogenase is highly sensitive to O₂. As such, for example, the activity of [Fe]-hydrogenase in C. reinhardtii is very sensitive to O₂ during H₂-photoproduction under photosynthetic conditions. Ghirardi et al., (1997) App. Biochem. Biotech. 63-65:141-151. Oxygen inhibition of[Fe]-hydrogenases is a major drawback in the use of green alga for H₂ production.

One approach to overcoming this H₂ production limitation is to stress the C. reinhardtii under photoheterotrophic, sulfur-deprived conditions that minimize O₂-photoproduction levels and result in sustained H₂-production. However, this approach does not result in optimal yields and requires the use of suilir-deprived/oxygen limited production techniques. Recently, CO and O₂ inhibition of hydrogenase activity in alga has been focused on the putative role of the H₂-channel. For example, it has been shown that the positioning of the Fe₂-atom in the enzyme's active site is directly at the active-site/H₂-channel interphase, where it is easily accessed by either CO or O₂ diffusing through the channel. Lemon et al., (1999) Biochem. 38:12969-12973; Bennett et al., (2000) Biochem. 39:7455-7460. Further, a naturally occurring O₂-resistant [NiFe]-hydrogenase has been shown to have a narrower active site/H₂-channel interphase than the naturally occurring hydrogenase counterpart. Volbeda et al. (2002), Supra.

Against this backdrop the present invention has been developed.

SUMMARY OF THE INVENTION

The present invention provides oxygen-resistant hydrogenases for use in the bulk production of H₂ in green algae cultures. In a preferred embodiment, homology modeling between known hydrogenases, eg., CpI, and target hydrogenases, e.g., HydA1, was used to design and in silico engineer an oxygen-resistant [Fe]-hydrogenase having a reduced diameter H₂-channel. Constructed polynucleotides that encode oxygen-resistant [Fe]-hydrogenase enzymes are used to transform target host cells which were subsequently used in the photoproduction of H₂. In preferred embodiments, the target host cells are C. reinhardtii. The invention provides a solution to the problem of H₂ production by green algae when O₂ is present in the environment.

The present invention also provides host cells expressing oxygen-resistant [Fe]-hydrogenase. Host cells expressing the oxygen-resistant [Fe]-hydrogenase have significantly increased H₂ production, in the presence of O₂, as compared to similarly treated cells that do not express oxygen-resistant [Fe]-hydrogenase.

The present invention also provides polynucleotide molecules encoding HydA1V240W and other like oxygen-resistant hydrogenase polypeptides. The invention includes nucleic acid molecules that hybridize under high stringency conditions to the HydA1V240W polynucleotides (and other like oxygen-resistant hydrogenase polynucleotides) of the present invention. The invention also includes variants and derivatives of the oxygen-resistant [Fe]-hydrogeanse polypeptides, including fusion proteins that confer a desired function. The invention also provides vectors, plasmids, expression systems, host cells and the like, containing the oxygen-resistant [Fe]-hydrogenase of the invention.

These and various other features and advantages of the invention will be apparent from a reading of the following detailed description and a review of the appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1A shows a protein alignment using PILEUP/GENEDOC program. Amino acid residues highlighted in black represent identities between at least 5 of the iron-hydrogenases, and those highlighted in grey show similarity between at least 5 of the sequences (1A). FIG. 1B shows a theoretical structure of HydA1 using homology modeling to the solved X-ray structure of CpI. The left panel shows an overlay of HydA1 and CpI, with locations of the H₂-channels and the active sites, while the right panel shows the HydA1 structure.

FIG. 2 shows the protein sequence of HydA1 aligned to the catalytic core region of CpI. The sequences that form the H₂-channel domain are shaded either gray (similar) or black (identical).

FIG. 3A and 3B show how mutations made to the HydA1 H₂-channel result in predicted H₂-channel structures. Wild type HydA1 H2-channel structure is shown in 3B, while mutant H2-channel structures are shown in 3A. Note that for reference purposes, the channel has been divided into four zones (black line numbered 1-4).

FIG. 4A and 4B show side-orientation views from the active site (left) to the protein surface (right) of the H₂-channel of the wild type HydA1 (4A) and mutant HydA1V240W (4B).

FIG. 5 illustrates PCR products from C. reinhardtii HydA1V240W transformants mt18 and mt28. Genomic DNA isolated from cc849 (wild type), mt18, and mt28 were digested with either SacI (lanes 1-3) or EcoRI (lanes 4-6) and used as template in a PCR reaction with HydA1 specific primers. Lanes 1 and 4 are wild type, lanes 2 and 5 are mt18 and lanes 3 and 6 are mt28. Note that lane 7 is a DNA size marker and lane 8 is a pAExBle control. The upper band in the stained agarose gel corresponds to the HydA1 genomic copy, and the lower band corresponds to the HydA1 cDNA insert.

FIG. 6 illustrates hydrogenase activity as measured by the rate of H2 evolved (μmol H₂/mg ch1⁻¹/h⁻¹) under variant O₂ concentrations (0 to 3.5% final O₂ concentration) and plotted relative to the activity value obtained under completely anaerobic conditions.

FIG. 7 shows the activity of O₂-resistant [Fe]-hydrogenase as measured in a reduced MV assay. Samples of induced cells were taken and assayed for hydrogenase activity following exposure to various levels of O₂ (0-4% final O₂ concentration). Note that hydrogenase activity was measured as the rate of H2 evolved (μmol H₂/mg ch1⁻¹/h⁻¹) over a 30-minute incubation period and plotted relative to the activity value obtained under completely anaerobic conditions.

FIG. 8 shows a plasmid map for the plasmid pLam91-1.

FIG. 9 shows a plasmid map for the plasmid pA1ExBle.

DETAILED DESCRIPTION OF THE INVENTION DEFINITIONS

The following definitions are provided to facilitate understanding of certain terms

“Amino acid” or “residues” refers to any of the twenty naturally occurring amino acids as well as any modified amino acid sequences. Modifications may include natural processes such as posttranslational processing, or may include chemical modifications which are known in the art. Modifications include but are not limited to: phosphorylation, ubiquitination, acetylation, amidation, glycosylation, covalent attachment of flavin, ADP-ribosylation, cross-linking, iodination, methylation, and alike. Amino acid residue characterization can be found in numerous citations, for example Stryer, 1995, Biochemistry, throughout the text and 17-44.

“Expression” refers to transcription and translation occurring within a host cell. The level of expression of a DNA molecule in a host cell may be determined on the basis of either the amount or corresponding mRNA that is present within the cell or the amount of DNA molecule encoded protein produced by the host cell (Sambrook et al., 1989, Molecular Cloning: A Laboratory Manual, 18.1-18.88).

“Genetically engineered” refers to any recombinant DNA or RNA method used to create a host cell that expresses a target protein at elevated levels, at lowered levels, or in a mutated form. Typically, the host cell has been transfected, transformed, or transduced with a recombinant polynucleotide molecule, and thereby been altered so as to cause the cell to alter expression of the desired protein. Methods for genetically engineering host cells are well known in the art. (See Current Protocols in Molecular Biology, Ausubel et al., eds. (Wiley & Sons, New York, 1988, and quarterly updates)). Genetically engineering techniques include but are not limited to expression vectors, targeted homologous recombination and gene activation (see, for example U.S. Pat. No. 5,272,071 to Chappel) and trans activation by engineered transcription factors (See Segal et al., 1999, Proc Natl Acad Sci USA 96(6):2758-63).

“Hybridization” refers to the pairing of complementary polynucleotides during an annealing process. The strength of hybridization between two polynucleotide molecules is impacted by the homology between the two molecules, stringency of the conditions involved, and melting temperatures of the formed hybrid and the G:C ratio within the polynucleotide. For purposes of the present invention stringency hybridization conditions refers to the temperature, ionic strength, solvents, etc, under which hybridization between polynucleotides occurs.

“Identity” refers to a comparison between pairs of nucleic acid or amino acid molecules. Methods for determining sequence identity are known in the art. For example, computer programs have been developed to perform the comparison, such as the Gap program (Wisconsin Sequence Analysis Package, Version 8 for Unix, Genetics Computer Group, University Research Park, Madison Wis.), that uses the algorithm of Smith and Waterman (1981) Adv Appl Math 2:482-489.

“Isolating” refers to a process for separating a nucleic acid or polypeptides from at least one contaminant with which it is normally associated. In preferred embodiments, isolating refers to separating a nucleic acid or polypeptide from at least 50% of the contaminants with which it is normally associated, and more preferably from at least 75% of the contaminants with which it is normally associated.

The term “nucleic acid” refers to a linear sequence of nucleotides. The nucleotides are either a linear sequence of polyribonucleotides or polydeoxyribonucleotides, or a mixture of both. Examples of nucleic acid in the context of the present invention include—single and double stranded DNA, single and double stranded RNA, and hybrid molecules that have mixtures of single and double stranded DNA and RNA. Further, the nucleic acids of the present invention may have one or more modified C) nucleotides.

The term “PCR” or “polymerase chain reaction” refers to the process to amplify nucleic acids as described in U.S. Patent Nos. 4,683,105 and 4,683,202, both owned by Roche Molecular.

“Host cell” refers to cells containing a target nucleic acid molecule, for example a heterologous nucleic acid molecule such as a plasmid or other low molecular weight nucleic acid, in which case the host cell is typically suitable for replicating the nucleic acid molecule of interest. Examples of suitable host cells useful in the present

E. Coli DH5α cells, as well as various other bacterial cell sources, for example the E. Coli strains: DH10b cells, XL1Blue cells, XL2Blue cells, Top10 cells, BB101 cells, and DH12S cells, yeast host cells from the genera including Saccharomyces, Pichia, and Kluveromyces and green alga, for example Chlamydomonas reinhardtii.

“Hybridization” refers to the pairing of complementary polynucleotides during an annealing period. The strength of hybridization between two polynucleotides molecules is impacted by the homology between the two molecules, stringent conditions involved, the melting temperature of the formed hybrid and the G:C ratio within the polynucleotides. High stringency conditions include, for example, 42° C., 6×SSC, 0.1% SDS for 2 hours.

“Nucleic acid” or “NA” refers to both a deoxyribonucleic acid and a ribonucleic acid. As used herein, “nucleic acid sequence” refers to the order or sequence of deoxyribonucleotides or ribonucleotides along a strand. They may be natural or artificial sequences, and in particular genomic DNA (gDNA), complementary DNA (cDNA), messenger RNA (mRNA), transfer RNA (tRNA), ribosomal RNA (rRNA), hybrid sequences or synthetic or semisynthetic sequences, oligonucleotides which are modified or otherwise. These nucleic acids may be of human, animal, plant, bacterial or viral origin and the like. They may be obtained by any technique known to persons skilled in the art, and in particular by the screening of libraries, by chemical synthesis or by mixed methods including the chemical or enzymatic modification of sequences obtained by the screening of libraries. They may be chemically modified, e.g. they may be pseudonucleic acids (PNA), oligonucleotides modified by various chemical bonds (for example phosphorothioate or methyl phosphonate), or alternatively oligonucleotides which are functionalized, e.g. which are coupled with one or more molecules having distinct characteristic properties. In the case of deoxyribonucleic acids, they may be single- or double-stranded, as well as short oligonucleotides or longer sequences. In particular, the nucleic acids advantageously consist of plasmids, vectors, episomes, expression cassettes and the like. These deoxyribonucleic acids may carry genes of therapeutic interest, sequences for regulating transcription or replication, anti-sense sequences which are modified or otherwise, regions for binding to other cellular components, and the like.

“Oxygen resistant” refers to any measurable decrease in oxygen sensitivity in a hydrogenase as compared to a hydrogenase having a reference oxygen sensitivity, for example, as compared to a wild type hydrogenase from which an oxygen-resistant hydrogenase enzyme has been made.

“Oxygen sensitive” refers to the wild type or reference oxygen sensitivity found in a native hydrogenase.

“Protein,” “peptide,” and “polypeptide” are used interchangeably to denote an amino acid polymer or a set of two or more interacting or bound amino acid polymers.

Green Algae and Iron Hydrogenase

Green algae, e.g., Chlamydomonas reinhardtii, cultured under anaerobic conditions synthesize an enzyme known as iron-hydrogenase ([Fe]-hydrogenase). As shown in FIG. 1A, several representative [Fe]-hydrogenase enzymes are aligned showing general sequence identity for the family of proteins. Generally, overall sequence identity for [Fe]-hydrogenase family members is usually at least 45% and sequence identity for the H₂-channel between family members is at least 66%.

In general, [Fe]-hydrogenase enzymes characteristically possess a catalytic site consisting of a bimetallic center containing two Fe atoms (2Fe-center), bridged by cysteinyl sulfur to an electron relay [4Fe4S] center (4Fe-center). The iron atoms of the catalytic 2Fe-center are joined together by a combination of organic, sulfur, and carbon monoxide ligands. The chemistry of the [Fe]-hydrogenase catalytic core is reactive with respect to hydrogen, typically possessing very high hydrogen-production rates. However, this same catalytic core is also highly sensitive to inactivation by oxygen. As a protective measure against inactivation by oxygen or other like molecules, the catalytic core is typically buried deep within the protein, where access to the core is limited. As a result, interface of the hydrogenase catalytic site with surface surroundings is principally limited to a single channel, termed the H₂-channel, that directs diffusion of synthesized hydrogen from the enzyme interior to the external environment. The H₂-channel is also the primary access route of oxygen to the metallo-catalytic site within hydrogenase enzyme. Reverse diffusion of the oxygen from the surface of the enzyme into the H₂-channel and on to the active site, allows oxygen to bind to the 2Fe-center, inactivating the enzyme. Under normal physiologic conditions this represents a fairly normal inhibitory response for the hydrogenase enzyme, however, under the artificial conditions of expressing bulk amounts of H₂, this is a fairly major limitation.

The present invention provides for the modification of the H₂-channel to reduce oxygen diffusion from the external environment to the enzyme's catalytic core. The present invention provides modifications to the H₂-channel that act as oxygen filters, preventing or reducing the diffusion of oxygen to the catalytic site within the hydrogenase enzyme. These modifications are at the same time insufficient to limit the ability of H₂ to diffuse out of the enzyme through the H₂-channel. Several mechanisms for the reduction of oxygen diffusion to the active site within the hydrogenase enzyme are provided, including targeted replacement of residues that line the H₂-channel with bulkier residues, so as to shield the 2Fe-center and/or reduce the diameter of the H₂-channel. In particular, the residues that line the H₂-channel are replaced with bulkier, hydrophobic residues, for example tryptophan or phenylalanine, so as to shield the 2Fe-center, as well as to reduce the size or volume of the catalytic site-H₂-channel interface. In addition, modifications to residues on the channel interior that approach and define the channel-solvent boundary (see portion 4 of FIG. 3).

In a preferred embodiment of the invention, a process for designing and engineering oxygen-resistant iron-hydrogenases has been developed. The engineering scheme targets the structure and or environment of the H₂-channel within the target hydrogenase, which is altered to be more selective in allowing the outward diffusion of hydrogen while simultaneously filtering our surface oxygen. Note that size-limited diffusion has been successfully used to generate filters for commercial use in the separation of gases, including the separation of hydrogen from oxygen. Menoff T. M.,

The present invention provides host cells for the expression of nucleic acid molecules for encoding an oxygen-resistant iron-hydrogenase, for example, C. reinhardtii that expresses an oxygen-resistant HydA1 or cyanobacteria that also expresses an oxygen-resistant HydA1. Example oxygen-resistant hydrogenases designed and engineered by the method of the present invention include V240W, A78W, A244W, A248W, G86W, and L93W. These oxygen-resistant hydrogenase enzymes have the same primary structure as HydA1 with the exception that A at residue 78 is replaced with W. Note that other residues besides W, including synthetic and derivatized amino acids, are envisioned for substitution into the H₂-channel, as long as they limit O₂ diffusion through the channel and allow H₂ diffusion out of the channel.

In addition, the invention provides for the bulk photoproduction of H₂ using the transformed host cells of the invention.

Identification of The Residues That Form the H₂-Channel Of A [Fe]-Hydrogenase

The iron-hydrogenase family of enzymes is a group of enzymes expressed in algae for metabolism of hydrogen. Iron-hydrogenase family members have been shown to have three distinct motifs that contain highly conserved residues, including a series of identifiable cysteine residues. Vignais et al., (2001) FEMS Micro. Rev. 25:455-501. In particular, motif 1 has the amino acid sequence PMFTSCCPxW, motif 2 has an amino acid sequence MPCxxKxxExxR and motif 3 has an amino acid sequence of FxExMACxGxCV. These three motifs have been identified in all iron-hydrogenase family members to date. The cysteine residues have been shown to either ligate the catalytic [4Fe-4S] center, or bridge the [4Fe-4S] to the [2Fe-2S] center, and there presence within the primary structure of the enzyme is highly conserved. One of the most studied iron-hydrogenase enzymes is CpI, having its primary, secondary and tertiary structures determined. Peters et al., (1998) Science 282:1853-1858. In preferred embodiments, CpI or other like known iron-hydrogenase enzymes, can be used in the design and engineering of oxygen-resistant hydrogenases (see below and FIG. 1A for potential iron-hydrogenase enzymes).

To identify the H₂-channel within a target hydrogenase, i.e., a polypeptide containing motifs 1-3 above, the primary sequence of the target hydrogenase must be compared to the primary sequence of a known hydrogenase. Once the two sequences have been aligned a level of identity is determined (see FIG. 1A and 2). Stothard P., (2000) BioTechniques 28(6) 1102 (hereby incorporated by reference in its entirety). For purposes of the present invention an overall identity of approximately 40%-45% or better should be found for the target hydrogenase. Further, an analysis of the target polypeptide's primary sequence is performed to predict the sequences that share homology with the H₂-chanel forming regions of other known iron-hydrogenases (similar patterns of residues that have been shown previously to form hydrophobic cavities). Montet et al., (1997) Nat. Struc. Biol; 4:523-526. It should be noted that because the H₂-channel is a conserved domain within all hydrogenases, other non-iron hydrogenase sequences can be used to identify the target hydrogenase H₂-channel. There should be at least 40%-45% identity between the known and unknown sequence between the H₂-channel sequences of the know and unknown hydrogenases. Once the region within the target polypeptide for the H₂-channel has been located, the channel is modeled into a three-dimensional structure showing the orientation of residues in relation to the channel and active site. Guex et al., (1997) Electrophoresis 18:27142723. (see below) In some embodiments, the analysis is extended to identify the residues corresponding to the active site within the target hydrogenase. Note that the active site of the target or unknown hydrogenase should share at least 90% homology for motifs 1-3, and in preferred embodiments shown complete identity with motifs 1-3 (see above). The combination of primary and tertiary structures of the target hydrogenase are compared to evaluate the identification of candidate regions for the final verification of the hydrogen-channel.

Methods For Designing and Engineering Oxygen-resistant Iron-Hydrogenases

As noted above, the present invention provides a model for generating a theoretical structure of a target H₂-channel within a target hydrogenase enzyme. In one embodiment, the theoretical structure is generated by homology modeling (see above) to the solved structure of other known [Fe]-hydrogenases, for example CpI. (see

hydrogenase active site and H₂-channel, and in other embodiments the homology modeling can be limited to the known hydrogenase H₂-channel sub-domains. A percent homology of the known hydrogenase (both identity and similarity) can be used to determined residue identity and similarity for the entire enzyme, the active site, the H₂-channel and the H₂-channel sub-domains (see overhead arrows in FIG. 2 and see discussion in previous section above). As such, the present invention provides a known hydrogenase based homology model that gives a reliable approximation of the target hydrogenase structure and H₂-channel environment. In a preferred embodiment, the known hydrogenase is CpI and the target hydrogenase is HydA1. Homology modeling can be performed using Swiss-model software as described in Guex et al. Guex et al (1997) Electrophoresis 18:2714-2723. Note, however, that other like programs can be used in this aspect, as is known in the art, e.g., Modeller program designed by Marti-Renom et al., (2000) Ann. Rev. Biophy, Biomol. Struct. 29:291-325; EsyPred3D designed by Lambert C. et al., (2002) Bioinformatics 18(9):1250-1256.

Typically, the homology modeling identifies the residues that project into the H₂-channel interior of the target hydrogenase. The channel environment is often composed of smaller hydrophobic residues, e.g., glycine, alanine, valine, but can contain phenylalanine and other like residues. For example, the H₂-channel of HydA1 contains mostly small hydrophobic residues with the exception of the larger phenylalanines at positions 252 and 355 (see FIG. 2, black dotted residues). A secondary structure is determined from the active site to the enzyme surface using the modeled structure above, and distances between side chain atoms of identified residues opposed to each other are determined. Guex et al., Supra An approximate average diameter of the channel over the distance from the catalytic site (Fe2-atom to the H-cluster [2Fe-2S]-center) to the protein surface is determined (see FIG. 3, 1-4) (typically by using the distances between the side chain atoms of opposed residues within the channel). It silico mutagenesis is performed on the identified hydrogenase H₂-channel structure to identify possible residues that can be modified to reduce the H₂-channel diameter. Mutagenesis criteria preferably involve conservative mutation of specific residues, selection of the lowest energy rotomer and energy minimization of the resulting structure using GROMOS. van Gunsteren, W. F. et al., (1996) Biomolecular Simulation, The GROMOS96 Manual and User Guide. Vdf Hochschulverlag ETHZ. Once an energy minimized structure is obtained, the dimensions of the target in silico mutagenized hydrogenase channel is determined. In preferred embodiments, one or more locations along the H₂-channel is designed via conservative mutation to be smaller in diameter than a corresponding non-mutated H₂-channel, typically this reduction is to a channel size of between approximately 5.0 and 2.4 Å in diameter, and preferably between 3.5 and 2.5 Å, a diameter that either limits or eliminates the ability of oxygen to diffuse through the modified H₂-channel. Note that the H₂-channel is in constant flux, as such diameter measurements are averages and not meant to be held to a static standard. Note that in embodiments of the present invention, more than one residue can be in silico mutated to design an optimum oxygen-resistant hydrogenase.

In an alternative embodiment, design of oxygen-resistant hydrogenase enzymes is provided by determining what substitutions/modifications of residues within the identified H₂-channel of a target [Fe]-hydrogenase can be performed to reduce the volume of the H₂-channel. Volume considerations include a reduction in the flow of gasses, i.e., O₂, through the channel in accordance with Stokes Einstein Equation and Fick's law.

Designed oxygen-resistant hydrogenases, having a reduced diameter H₂-channel, are genetically engineered and transformed into target host cells, for example, into C. reinhardtii, and tested for hydrogenase activity in the presence of O₂ via a modified Clark electrode or other known assay(s). In preferred embodiments, the oxygen-resistant hydrogenase is generated via site-directed mutagenesis. For example, to generate HydA1 mutants, the HydA1 gene of pA1ExBle can be mutagenized in vitro using the Quick Change XL Site-Directed Mutagenesis Kit (Stratagene). Host cells that have incorporated the designed enzymes having reduced oxygen sensitivity) can be used to photoproduce H₂ in an oxygen containing environment. Note that these host cells can also be treated with mRNA interference to repress the expression of native hydrogenases, while continuing to allow expression of the inventive engineered

Steered Molecular Dynamics (SMD)

In one embodiment, the in silico designed oxygen-resistant hydrogenase enzymes can be further analyzed for changed or reduced oxygen diffusion within their H₂-channel by applying SMD via the NAMD program. Kale L. et al., (1999) Computational Physics 151:283; Isralewitz B., (2001) Curr. Opin. Struc. Biol. 11:224. SMD analysis, therefore, provides confirmation and additional baseline data as to the efficiency of the channel modifications and their effects on O₂ diffusion within the proposed oxygen-resistant hydrogenase.

Oxygen-Resistant Hydrogenase Polypeptides

Oxygen-resistant hydrogenase enzymes of the invention include all proteins that can be constructed from the in silico mutagenesis methods discussed above. For example, any polypeptide having a predicted reduction in hydrogen-channel diameter or volume, as determined by the methods of the invention, is envisioned to be within the scope of the present invention.

In addition, oxygen-resistant hydrogenase enzymes of the invention include isolated polypeptides having an amino acid sequence as shown in FIG. 2 (Cr HydA1), and having one or more substitutions at residues V240, A78, A244, A248, G86, and L93 (note that substitution by tryptophan and other like amino acids is envisioned, including synthetic or derivatized amino acids) (also included are substitutions shown in Tables 1 and 2). The invention includes variants and derivatives of these oxygen-resistant [Fe]-hydrogenase enzymes, including fragments, having substantial identity to these amino acid sequences, and that retain both hydrogenase activity and enhanced tolerance to oxygen (see Example 3 for assays to determine hydrogenase activity in the presence of oxygen). In a preferred embodiment, the oxygen-resistant hydrogenase enzyme is HydA1V240W. Derivatives of the oxygen-resistant hydrogenases include, for example, oxygen-resistant HydA1 enzymes modified by covalent or aggregative conjugation with other chemical moieties, such as lipids, acetyl groups, glycosyl groups, and the like.

Oxygen-resistant hydrogenase enzymes of the present invention can be fused to heterologous polypeptides to facilitate purification. Many available heterologous peptides allow selective binding of the fusion protein to a binding partner, for example, 6-His, thioredoxin, hemaglutinin, GST, and the like.

Polypeptide fragments of the modified oxygen-resistant hydrogenase H₂-channel polypeptide (that include the relevant residue modification) can be used to generate specific anti-oxygen-resistant hydrogenase antibodies (monoclonal or polyclonal). Generated antibodies can be used to selectively identify expression of oxygen-resistant hydrogenases or in other known molecular and/or biochemical techniques, for example, in immunoprecipitation or Western blotting.

Variant oxygen-resistant hydrogenase enzymes include fusion proteins formed of a oxygen-resistant hydrogenase and a heterologous polypeptide. Preferred heterologous polypeptides include those that facilitate purification, stability or secretion.

Oxygen-Resistant Hydrogena se Polynucleotides, Vectors and Host Cells

The invention also provides polynucleotide molecules encoding the oxygen-resistant polypeptides of the invention. The polynucleotide molecules of the invention can be cDNA, chemically synthesized DNA, DNA amplified by PCR, RNA or combinations thereof.

The present invention also provides vectors containing the polynucleotide molecules of the invention, as well as host cells transformed with such vectors. Any of the polynucleotide molecules of the invention may be contained in a vector, which generally includes a selectable marker, and an origin of replication, for propogation in a host. The vectors also include suitable transcriptional or translational regulatory sequences, such as those derived from algae operably linked to the oxygen-resistant hydrogenase polynucleotide molecule. Examples of such regulatory sequences include transcriptional promoters, operators, enhances, and mRNA binding sites. Nucleotide sequences are operably linked when the regulatory sequence functionally relates to the DNA encoding the target protein. Thus, a promoter nucleotide sequence is operably linked to a oxygen-resistant hydrogenase DNA sequence if the promoter nucleotide sequence directs the transcription of the oxygen-resistant hydrogenase sequence.

Selection of suitable vectors for the cloning of oxygen-resistant hydrogenase polynucleotides of the invention will depend on the host cell in which the vector will be transformed/expressed. For example, the plasmid pLam91-1 (see FIG. 8) was used to create a 980 bp HydA1 PstI promoter fragment cloned into a unique PstI site of pLam91-1, creating the HydA1 promoter-HydA1 cDNA fusion construct, pA1Ax. The Ble^(r) cassette of pSP108 confers Bleomycin resistance in transformed C. reinhardtii, and was inserted into the Tfi 1 site of pA1Ax, creating pA1AxBle (see FIG. 9). This was particularly useful in the construction of expression oxygen-resistant [Fe]-hydrogenase vectors for use in green algae.

Suitable host cells for expression of target polypeptides of the invention include green algae, for example C. reinhardtii cells and cyanobacteria, both of which utilize water in growth, which is also a substrate for the hydrogenase enzymes. Typically, green algae cells are transformed by a glass bead method as is known in the art. Cells exhibiting the target selectable marker, for example resistance to bleomycin, are picked and patched onto fresh TAP+Ble plates and re-patched an additional 2-3 times to ensure the isolation of stable integrates.

H₂ Production

Green algal cultures that express oxygen-resistant hydrogenase of the invention may be used to photoproduce H₂ in the presence of oxygen. In one embodiment of the invention, the transformed cells are grown in a photobioreactor photoautotrophically, photoheterotrophically in TAP, or other like growth media to a concentration of 5-50 μg/ml chlorophyll, and H₂ harvested. Note that in some embodiments, the cells are grown under selective pressure that ensures that the cells maintain the oxygen-resistant hydrogenase, for example in bleomycin, where the construct used to transform the host cell confers the selective pressure.

In another embodiment, the oxygen-resistant hydrogenase of the invention may be transformed into target algae, under the control of the endogenous HydA1 promoter, for nighttime enzyme generation and daytime H₂-production. See Boichenko et al., (2003) Photoconversion of Solar Energy, Molecular to Global Photosynthesis: In Press.

It is envisioned that the proceeding discussion on the design, engineering, and construction of oxygen-resistant hydrogenases, as well as the subsequent tansformation of host cells with the designed hydrogenases, can be expanded to any iron hydrogenase known or identified in the future having the characteristics for iron hydrogenase enzymes discussed herein.

Having generally described the invention, the same will be more readily understood by reference to the following examples, which are provided by way of illustration and are not intended as limiting.

EXAMPLES Example 1 Computer Modeling of Hydrogenase H₂-Channel For Design of Oxygen-Resistant Hydrogenase Enzymes

To facilitate the design and engineering of mutant oxygen-resistant HydA1 enzymes, a theoretical structure of HydA1 was generated by homology modeling to the solved X-ray structure of Clostridium pasteuraianuin [Fe]-hydrogenase, CpI (FIG. 1B). The theoretical model was generated by homology modeling using Swiss-model software as described by Guex et al. Guex et al., (1997) Electrophoresis 18:2714-2723. The resulting HydA1 model was subjected to several rounds of energy minimization using GROMOS. An alignment of the HydA1 and CpI amino acid sequences show they share a high degree of homology (45% identity, 58% similarity) within the essential domains, i.e., active site and H₂-channel, that comprise the core region of [Fe]-hydrogenases (see FIG. 2). Stothard P., (2000) BioTechniques 28(6) 1102. Note that the degree of conservation increases for H₂-channel sub-domains, where the two proteins share 62% identity and 92% similarity (FIG. 2, overhead

HydA1 homology model provides a reasonable approximation of the HydA1 structure and the H₂-channel environment.

The detailed study of the HydA1 H₂-channel structure was performed, at least partly, to identify residues that project into the H₂-channel interior. In general, the channel environment was primarily composed of smaller hydrophobic residues, e.g., glycine, alanine, valine, with the exception of the larger phenylalanines at positions 252 and 355 (FIG. 2, black dotted residues). The secondary structure of the H₂-channel was organized into two α-helices and two β-sheets, which extend from the active site to the enzyme surface. The distance between side chain atoms of residues that oppose each other were measured to approximate the average channel diameter over the distance from the catalytic site (Fe2-atom of the H-cluster [2Fe-2S]-center) to the protein surface (1 to 4, FIG. 3). The channel measured 3.85 to 7.44 Å in diameter over a distance of 24 to 27 Å, making the channel diameter greater than the effective diameters of both H₂ (2.8 Å) and O₂ (3.5 Å). As a result, the predicted size of the HydA1 H₂-channel is sufficient to function in H2 diffusion from the active site to the surface, but it is also large enough to allow for the inward diffusion of the inhibitor O₂.

These results suggest that engineering O₂ tolerance into HydA1 might be accomplished by altering the residues that line the interior of the channel so as to reduce the diameter of the channel and thereby limit O₂ diffusion to the active site. The potential to reduce the diameter of the channel via residue substitution was initially tested in silico by mutating the H₂-channel of the HydA1 model. The mutagenesis criteria involved conservative mutation, i.e., hydrophobic→hydrophobic, of specific residues, selection of the lowest energy rotomer, and energy minimization of the resulting structure using GROMOS. Once an energy-minimized structure was obtained, the dimensions of its channel was determined. Several of the channel residues proved to be unameanable to mutation and were left unchanged, i.e., 182, L89, F252 and F355, i.e., the Guex program determined that changes at these locations would provide only minimal (non-significant) change to the H₂-channel diameter/volume. However, promising mutants were generated from alteration of several residues that were spaced over the entire length of the channel (see FIG. 2 and Table 1). Substantial reductions of channel diameter were obtained by mutating residues A78 and V240 (proximal to active site); AM44, A248 and G86(mid-channel);

and L93 (protein-solvent boundary, distal to active site) to bulkier amino acids (Table 1 and FIG. 3). The individual mutations listed in Table 1 caused reductions in diameter that ranged from 0.5 to 1.90 Å (Table 2). The HydA1 mutant that combined the A248I and L93F mutations located at the channel-solvent boundary (FIG. 3, zone 4) showed an average decrease in size from 5.21 to 3.34 ÅA, less than the effective diameter of O₂ (3.5 Å). When the mutations listed in Table 2 were combined into a single HydA1 mutant, the average overall channel diameter was reduced from an average 5.71 to an average of 4.31 Å (Table 2), noting however that there are several locations along the H₂-channel with reductions in the diameter at or near the average diameter of O₂. TABLE 1 Predicted Effects of Selected HydA1 H2-Channel Mutations on Channel Environment Mutation Location Effects A781 Adjacent to Fe₂-atom, Bulkier isoleucine side across from V240 chain projects closer to V240, and the [2Fe-2S]- center Fe₂ atom. V240W Adjacent to A78, above Bulkier tryptophan side Fe₂-atom. chain reduces distance to A78, and partially shields Fe₂-atom. A244L Mid-channel, opposes I82 Bulkier leucine side chain projects further into channel. G86I Mid-channel, opposes Isoleucine side chain adds A248 bulk, and projects into channel. A248I Mid-channel, near surface, Isoleucine extends further opposes G86 and L89 into channel, adds more bulk to hydrophobic surface. L93F Channel-Surface boundary Narrows the channel opening at protein surface- solvent boundary

TABLE 2 Distances Between Channel Determinants In HydA1 and HydA1 Mutants Based on Modeling Studies Distances (Å) Average Zone Size (Å) Determinant HydA1 HydA1 HydA1 HydA1 Zone^(a) Pairs^(b) wild type mutant^(c) wild type mutant^(c) 1 A(I)78::Fe₂ 6.25 4.66 6.90 5.04 V(W)240::Fe₂ 6.14 5.02 V(W)240::A(I)78 7.43 5.43 2 A(L)244::I82 4.50 3.90 5.30 4.87 F355::G(I)86 6.86 5.80 3 G(I)86::A(I)248 4.38 3.54 5.42 3.92 G(I)86::T247 6.01 4.39 A(I)248::L89 3.85 3.26 A(I)248::F355 7.44 4.54 4 L90::A(I)248 6.11 3.67 5.21 3.34 L(F)93::F252 4.31 3.01 ^(a)The locations of H2-channel zones are identified in FIG. 3. ^(b)Determinants are identified as wild-type, with corresponding mutations in parentheses. ^(c)Measurements are averages of a HydA1 mutant possessing all the identified mutations within the designated zone.

The above results indicate that modeling of the HydA1 structure has revealed a hydrophobic channel extending from the active site to the enzyme surface. This channel would appear to be conserved in other [Fe]-hydrogenases. The channel's secondary structure is mainly a-helical, which suggests that the channel domain is fairly rigid. Perhaps the rigidity of the channel structure helps to prevent its collapse during folding. Volbeda et al., (2002) Int. J. Hyd. Energy 27:1449-1461. Rigidity would also be expected to contribute to conformational stability of the channel in the folded protein, and a static model should give reasonable approximations of shape and size. Our measurements of the HydA1 channel demonstrate that it is sufficient in diameter not only to allow for diffusion of the product H₂ but also the larger-sized inhibitors O₂ and CO. Since enzyme inhibition occurs quickly (minutes), following exposure of O₂ (Happe et al., (1994) Eur. J. Biochem. 222:769-774), the channel with our analysis. This data illustrates the utility of the present invention for engineering O₂-resistant, [Fe]-hydrogenase by manipulation of residues within the conserved H₂ channel. This modeling approach can be used in enzymes that have one channel or multiple channels to reduce inhibitor access to an enzyme active site.

Example 2 C. reinhardtii Can Be Transformed with HydA1 H₂-Channel Mutants

To test the ability of the predicted HydA1 H₂-channel mutants for limiting O₂ inhibition, an algal HydA1 expression system was created using the HydA1 endogenous promoter. From the modeling discussed in Example 1, the V240W mutation was selected for further examination. In vivo expression of the V240W mutant was performed and further testing of the mutant for O₂ resistance hydrogenase activity performed. Note that the V240W mutation is predicted to cause a constriction of the channel near the active site (see FIG. 4). In addition, the tryptophan projects over the Fe₂-atom, partially shielding it from the channel domain.

The Chilamydomonas reinhardtii strain cc849 (cw10, mt-) was used as the wild type parent strain throughout the remainder of this Example. Growth of liquid cultures were performed photoheterotrophically in TAP medium (Harris E, (1989) The Chlamydomonas Source Book, Academic Press, New York) with a continuous stream of 5% CO₂ under cool-white fluorescent light (150 μE/m⁻²/s⁻¹ PAR). Growth on solid medium was performed on TAP agar plates (TAP medium with 1.4% w/v agar). Note that when selection of Bleomycin resistance was performed, solid TAP medium was supplemented with 10 μg/ml Zeocin (Invitrogen).

A plasmid construct pLam91-1, containing the HydA1 cDNA and 3′-terminator regions cloned into the EcoRI-XhoI sites of pBluescript SK, was used to generate an algal HydA1 expression construct. A 980 bp HydA1 PstI promoter fragment was cloned into the unique PstI site of Lam91-1, creating the HydA1 promoter-HydA1 cDNA fusion construct, pA1Ex. The Ble^(r) cassette of pSP108 that confers Bleomycin resistance in transformed C. reinhardtii (Stevens et al., (1996) Mol Gen Genet 251:23-30) was inserted into the TfiI site of pA1Ex, creating pA1ExBle.

Site-directed mutagenesis was performed on HydA1 to generate HydA1 mutants for expression in C. reinhardtii. The HydA1 gene pA1ExBle was mutagenized in vitro using the Quick Change XL Site-Directed Mutagenesis Kit of Stratagene. Oligonucleotides (Integrated DNA Technologies) used for mutagenesis were designed based on the kit requirements. Mutant pA1ExBle constructs were sequenced to confirm the presence of individual mutations. The HydA1 mutant, V240W, contains a valine to tryptophan substitution at amino acid position 240 of the mature protein.

C. reinhardtii cells were next transformed by the glass bead method as is known in the art (see also Harris E) using 10 μg of linearized pA1ExBleV240W DNA. Following transformation, cells were cultured overnight in 2 ml of TAP medium to allow for cell recovery and phenotypic expression of Ble^(r). Transformed cells were harvested by centrifugation (2000×g, 5 minutes), resuspended in 1.5 ml TAP soft agar (TAP with 0.8% w/v agar) and spread onto TAP+Ble agar plates. Plates were incubated in the light for a period of 1-2 weeks and Ble^(r) colonies picked. Resistant colonies were patched onto fresh TAP+Ble plates, and re-patched an additional 2-3 times to ensure the isolation of stable integrates.

To ensure that the HydA1 cDNA genomic insert having the V240W mutation was present in the transformed C. reinhardtii, PCR and sequencing was performed on Ble^(r) transformants. Total genomic DNA was isolated from individual Ble^(r) transformants using the Plant Genomic Kit (Qiagen). A total of 0.5 to 1.0 μg of purified genomic DNA was digested with either SacI or EcoRI and used as template in a PCR reaction consisting of the HydA1 internal primers (5′-CACGCTGTTTGGCATCGACCTGACCATCATG-3′ and 5′-GCCACGGCCACGCGGAATGTGATGCCGCCCC-3′), 1 unit KOD HotStart polymerase (Novagen), 10 mM MgSO4, 25 mM of each dNTP, 2% DMSO (v/v), and water to a total volume of 50 μl. The presence of a HydA1 cDNA genomic insert results in an additional 780 bp HydA1 cDNA product together with the 1120 bp HydA1 genomic product. PCR reactions were run on 1× TAE agarose gels (1.25% agarose w/v), stained with ethidium bromide, and photographed (not shown). The 780 bp band, corresponding to the HydA1 cDNA insert, was purified and sequenced to confirm the presence of V240W mutation. Two Ble^(r) C. reinhardtii clones, mt18 and mt28, were shown to possess the HydA1V240W construct (see FIG. 5).

Example 3 Green Alga, C. reinhardtii, Transformed with Oxygen-Resistant Hydrogenases Are Effective In The Bulk Production of H₂

The O₂-sensitivity of [Fe]-hydrogenase activity in strains mt18 and mt28 carrying the HydA1V240W mutation was tested in either whole cells or whole cell extracts of anaerobically induced cultures. Hydrogenase activities were measured as H₂ gas photoproduction by whole cells as previously described. Ghirardi et al., (1997) App. Biochem. Biotech. 63-65:141-151; Flynn et al., (2002) Int. J. Hyd Energy 27:1421-1430. Briefly, cells were grown photoheterotrophically in TAP to a concentration of 15-20 μg/ml chlorophyll, harvested and resuspended at 200 μg/ml chlorophyll in phosphate induction buffer. Ghirardi et al. Clark electrode measurement of O₂-resistant hydrogenase activity was performed by adjusting the O₂ concentration in the electrode chamber to a set level between 0% and 4%. Once the O₂ level had stabilized, a stream of Ar gas was passed over the chamber to maintain a constant O₂ concentration. A 0.2 ml sample of induced cell suspension was injected into the chamber, and the cells kept in the dark for a two minute period. Light dependent H₂-photoproduction activity was then induced by illumination.

In addition, to measure hydrogenase activity directly, reduced methyl viologen (MV) was used as an artificial electron donor for H₂ production by solubilizing whole cells as previously described. Flynn et al. Tolerance to O₂ was measured by incubating 1 ml of induced cells in a dark, sealed glass bottle and injecting O₂ to achieve a final atmosphere of 1 to 4% (v/v). Samples were incubated for two minutes then purged with Ar gas for five minutes. A 1 ml mixture of reduced MV and Triton X-100 in a phosphate buffer was added, samples were mixed for three to five minutes, and 0.1 ml of 100 mM reduced Na-dithionite injected to start the reaction. The reaction mixtures were incubated for 30 minutes at room temperature with stirring, and reactions were stopped by the addition of 0.1 ml 20% trichloroacetic acid (TCA). The hydrogen content of a 0.2 ml headspace sample was measured by gas chromatograph. Three separate headspace samples were assayed, and the values were averaged to attain final hydrogen-production rates.

As shown in Table 3, all three strains, cc849, mt18, and mt28, exhibited similar levels of hydrogenase activity (rate of H₂ photoproduction) under completely anaerobic conditions. Note that as has been shown in previous studies (Ghirardi et al, supra; Flynn et al, supra), pretreatment of induced wild type cells with O₂ is sufficient to cause a significant decline in H2 production rate (FIG. 6, white bars). When induced wild-type cells were pre-treated with O₂ at a concentration of 1.7 to 3.5%, the H2 photoproduction rate declined by 90 to 100% respectively. However, the exposure of mt18 or mt28 induced cells to similar O₂ treatments showed H2 photoproduction activity had significant resistance to inactivation. After exposure to 1.7 to 2.2% O₂ concentrations, the H₂ photoproduction rates remained 3.8 to 7 fold higher in mt18, and 3.2 to 13 fold higher in mt28 compared to activities in wild-type cells under identical conditions (see FIG. 6). At 3.5% O₂ treatment, the H2 photoproduction rates in both mt18 and mt 28 were low, but detectable, whereas residual activity in wild-type cells was undetectable (FIG. 6). TABLE 3 Hydrogenase Activity By the Clark Electrode Assay H₂ Photoproduction Rate Strain (μmol H2/mg chl⁻¹/h⁻¹) cc849 10.4 mt18 14.1 mt28 10.7

The light-induced production of hydrogen by whole cells is a metabolic process and depends on many electron transfer steps. Zhang et al., (2000) Trends Biotech. 18(12):506-511; Melis et al., (2001) Plant Physiol. 127:740-748; Melis et al., (2000) Plant Physiol. 122:127-135. A more direct measurement of hydrogenase activity can be accomplished in solubilized whole cells using reduced MV (Mv_(red)) as electron donor for H2 gas production by hydrogenase in the dark. Under completely anaerobic conditions, the Mv_(red)→H₂ reaction rates were similar in value for either induced wild-type or mutant cells (see Table 4). As shown in FIG. 7, a two-minute exposure of induced wild-type to various O₂ concentrations caused hydrogenase activity to decline. After exposure of O₂ concentrations of 1% to 4% hydrogenase activities in wild-type cells decreased to between 10 and 1.5% respectively (FIG. 7), similar to the results shown in FIG. 6. In comparison, both mtl6 and mt28 containing the HydA1V240W construct exhibited significant levels of O₂ resistant hydrogenase activity (see FIG. 7). Exposure of mt18 to O₂ at 1% to 4% concentration caused hydrogenase activities to decline by 76% to 96%, whereas mt28 activities declined only 12% to 76% (FIG. 7). As a result, mt18 hydrogenase activities were 2- to 3-fold higher, and mt28 activities 8- to 15-fold higher than activities in wild-type cells after exposure to similar O₂ treatments. TABLE 4 Hydrogenase Activity By the Methyl Viologen Assay H2 Photoproduction Rate Strain (μmol H2/mg chl⁻¹/h⁻¹) cc849 31.3 mt18 32.9 mt28 35.5

This Example illustrated the utility of modeling residue substitutions within the H₂-channel to constrict the channel from O₂ passage to the [2Fe-2S]-center. In particular, the Example illustrated that substitution of tryptophan for valine at position 240 of HydA1 caused an increase tolerance to O₂ in the mutant hydrogenase. The difference in the structure change made to HydA1V240W and the effects of that change are similar to the observed differences in structure and O₂-resistance of H₂-sensing [NiFe]-hydrogenases compared to catalytic [NiFe]-hydrogenases. Volbeda et al., (2002) Int. J. Hyd. Energy 27:1449-1461; Bernhard et al., (2001) 276:15592-15597. Active-site proximal channel residues of O₂-resistant, H₂-sensing [NiFe]-hydrogenases contain the bulky, hydrophobic amino acids isoleucine and phenylalanine. Identical positions in the O₂-sensitive, catalytic [NiFe]-hydrogenases encode the smaller-sized residues valine and leucine respectively. The difference in amino acid composition is suggested to result in the shielding of the [NiFe]-cluster and constriction of the channel. Volbeda supra and Bernhard supra.

The invention has been described with reference to specific examples. These examples are not meant to limit the invention in any way. It is understood for purposes of this disclosure, that various changes and modifications may be made to the invention that are well within the scope of the invention. Numerous other changes may be made which will readily suggest themselves to those skilled in the art and which are encompassed in the spirit of the invention disclosed herein and as defined in the appended claims.

This specification contains numerous citations to patents and publications. Each is hereby incorporated by reference in their entirety for all purposes. 

1. An oxygen-resistant iron hydrogenase derived from an oxygen sensitive iron hydrogenase by the substitution of one or more amino acid residues within the hydrogen channel of the oxygen-sensitive iron hydrogenase.
 2. The oxygen-resistant iron hydrogenase of claim 1 wherein said substitution is made to one or more amino acid residues selected from the group consisting of residues 78, 240, 244, 86, 248, 247, 82, 89, 355, 93 and 252 of HydA1 iron hydrogenase.
 3. The oxygen-resistant iron hydrogenase of claim 1 wherein the substituted amino acid which is at a residue within said hydrogen channel has a side chain volume which is larger than the side chain volume of the amino acid at the same residue in said oxygen-sensitive iron hydrogenases.
 4. A nucleic acid encoding the oxygen-resistant iron hydrogenase of claim
 1. 5. A vector comprising the nucleic acid of claim
 4. 6. Host cell transformed with the vector of claim
 5. 7. The host cell of claim 6 wherein the host cell is green algae.
 8. A process for producing hydrogen in green algae comprising culturing green algae cells transformed with the vector of claim 5 under conditions wherein the nucleic acid encoding said oxygen-resistant iron hydrogenase is expressed.
 9. The process of claim 8 wherein said green algae is Chlanydomonas reinhardtii.
 10. A method for making a nucleic acid encoding an oxygen-resistant iron hydrogenase comprising: comparing the sequence of a first oxygen-sensitive iron hydrogenase having known amino acid residues that form a hydrogen channel in said first hydrogenase with a second iron hydrogenase to identify residues that form a hydrogen channel in said second hydrogenase, and forming a nucleic acid encoding a derivative of said iron hydrogenase wherein one or more amino acid residues within said hydrogen channel are modified to reduce the oxygen sensitivity of said second hydrogenase.
 11. The method of claim 10 wherein the crystal structure of said first hydrogenase is known and said comparing is by in silico homology modeling.
 12. The method of claim 10 wherein said first hydrogenase is CpI and said second bydrogenase is HydA1.
 13. A process for making a green algae capable of hydrogen production in the presence of oxygen comprising transforming a green algae with the vector of claim
 5. 14. A method of making an oxygen-resistant iron-hydrogenase comprising: determining the diameter of an H₂-channel defined by a set of amino acid residues in an oxygen-sensitive iron hydrogenase by measuring the distance between the side chains of one or more of the amino acid residues of said set to identity one or more diameter determining amino acid residues; and modifying one or more of said diameter determining residues in said oxygen-sensitive iron hydrogenase to reduce the effective diameter of said H₂-channel to form an oxygen-resistant iron hydrogenase, whereby the diffusion of oxygen within said modified H₂-channel in said oxygen-resistant iron hydrogenase is reduced as compared to the diffusion of oxygen in said H₂-channel of said oxygen-sensitive iron hydrogenase.
 15. The method of claim 14, wherein said H₂-channel residues are identified by homology modeling to a known x-ray structure of one or more known hydrogenases.
 16. The method of claim 15, wherein substitution at one or more of said diameter determining residues of said H₂-channel is made in the form of an in silico molecule wherein the most energentically favorable rotamer for the substituted side chain is used to minimize the energy of at least the H₂-channel of said molecule.
 17. The method of claim 16 wherein GROMOS is used to compute said minimized energy of said molecule.
 18. The method of claim 16 wherein the diameter of said H₂-channel in said energy minimized molecule is compared to the diameter of said H₂-channel prior to in silico modification as a basis for selection of one or more mutations to be formed in said oxygen-sensitive iron hydrogenase protein. 